Spectral Analysis of Correlated One-Dimensional Systems with Impurities 
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An averaging procedure is proposed to account for spectral features of correlated one-dimensional 
systems in the presence of non-magnetic impurities. The dynamical spin structure factor for a 
corresponding random ensemble of Heisenberg chain segments is calculated by exact numerical 
diagonalization. It is shown that a few-pole approximation is sufficient to describe the numerical 
results. A similar analysis is proposed for the discussion of experimental spectra, such as obtained 
by inelastic neutron scattering measurements on Zn-doped CuO chains. By examination of the 
disorder-induced pseudo-gap, the loss of spectral weight, and the discrete peak structures due to 
smallest-cluster contributions, the underlying impurity distribution function can be determined. 



It is known that electrons in one dimension are local- 
ized by an infinitesimal amount of disorder. Analo- 
gously, their spin degrees of freedom develop a spectral 
gap when impurities are placed into an infinite chain. 
The branch cut corresponding to the collective spinon 
continuum of the pure system splits into a discrete set 
of poles characterized by conformal towers of excita- 
tions, and the spacing between the poles in a given one- 
dimensional segment decreases with increasing size. In 
this paper it is proposed that an average over a distribu- 
tion of one-dimensional segments, defined by the regions 
between randomly placed impurities, should be taken 
when calculating physical quantities. At low energies, the 
finite system gaps of all contributing segments then add 
up to form a pseudo-gap. By analyzing the peak struc- 
ture, the overall spectral weight, and the pseudo-gap fea- 
tures of quasi-one-dimensional crystals in the presence of 
non-magnetic impurities, such as Zn-doped CuO chains 
and non-stochiometric a'-NaV2 05, |^ the respective 
distribution function of lengths of chain segments can 
thus be deduced, indicating the extent to which the in- 
troduction of defects has pushed the system into a meso- 
scopic regime. 

As a non-trivial paradigm for correlated one- 
dimensional systems, let us thus focus on the antifer- 
romagnetic spin-1/2 Heisenberg chain in the presence of 
non-magnetic impurities. Such a model is realized for 
example in Zn-doped SrCu02, where the antiferromag- 
netic superexchange between neighboring Cu^+ dx2_y2 
electrons is mediated by the filled O^^ p-orbitals. By 
substituting Zn^+ for Cu^+ , static vacancies are created, 
and the infinite chain is separated into segments of length 
I which follow a distribution function P{1). 

The dynamical structure factor S(q,a;) of the infinite 
antiferromagnetic spin-1/2 Heisenberg chain is described 
approximately by a two-spinon continuum of excitations 
which can be probed by inelastic neutron scattering 
experiments. |^ It is bounded from below by the des 
Cloiseaux-Pearson dispersion, LOi{q) = (7rJ/2)|sin(q)|, ||] 
and from above by the maximum energy of two unbound 
spinous W2(q) = 7rJ|sin(q/2)|. Q Thus at q=7r the two 
spinous have maximal phase space, and the dynamical 



structure factor diverges as S(7r,w) ~ (l/a;)ln(l/a;) at 
the lower boundary. |^,^ For finite systems with an even 
number of sites, I, the corresponding branch cut splits 
into a discrete set of l/A poles for /=4,8,12,... , and 
(?+2)/4 poles for /=2,6,10,... (see Fig. 1(a)). Chains 
with an odd number of sites do not contribute at q=7r. 

The finite-size scaling behavior of the pole positions, 
shown in Fig. 1(b), is well described by the towers of 
excitations found in conformal field theory (CFT), when 
logarithmic corrections due to marginally relevant opera- 
tors are taken into account. |p|-pT[| Thus the pole positions 
are of the form 



Lu,{l)^a,/l + /3J{l\n{l)), 
where a,; and (3i are treated as fit parameters. 
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FIG. 1. (a) Finite-size scaling of the q=7r dynamical 
structure factor of the antiferromagnetic spin-1/2 Heisenberg 
chain. In the thermodynamic limit {I = oo) the discrete set 
of poles merges into a single branch cut. The (5-peaks have 
been given a width of e=0.1J. (b) Finite-size scaling of the 
first three pole positions in (a), (c) Finite-size scaling of the 
corresponding pole weights. The solid lines in (b) and (c) are 
fits to the finite cluster data discussed in the text. 

In Fig. 1(b), the first three pole positions for finite-size 
clusters up to ^=24 (symbols) are shown along with the 
fits to Eq. 1 (solid lines). There is good agreement of ai 
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I e.g. af 
0.94af^^. 



7rV/2, 



with the predictions from CFT, 
while the fit yields ai = 4.614J ; 

The pole positions may be calculated directly from the 
Bethe Ansatz equations or approximately using CFT, 
whereas their respective matrix elements are accessible 
only to numerical diagonalizations of small clusters. In 
Fig. 1(c) it is shown that a fit of the finite cluster pole 
weights, Ai{l), to the form 



Ai{l) = Ui + 6iln(l + l) + Ci ln(l + ln(l + I)), 



(2) 



gives good results, and thus permits extrapolation to 
larger system sizes. The dynamical spin structure fac- 
tor at q—TT for a system of length / can then be written 
in the form 



(3) 



where the sum is taken over all poles. In practice, the 
(5-peaks are replaced by Lorentzians of width e, which 
throughout this paper will be taken as e = O.IJ. 

Let us now turn to the effect of vacancies entering an 
ideal infinite chain i.e. by doping the system with non- 
magnetic impurities. Even mobile carriers which can be 
introduced into the chain by out-of-plane substitutions 
(such as La^+ for Sr2+ in SrCu02 Q) tend to become 
localized due to small random potentials, and hence may 
take the same role as static vacancies. If the impuri- 
ties enter the chain in a completely random manner, the 
lengths of the chain segments thereby created follow a 
distribution function 



P{1) = pexp{-pl), 



(4) 



where p is the concentration of impurities per lattice site 
(0 < p < 1). [|l3[ Note that P{1) is properly normalized, 
and that its nP^ moment is given by (P) = n!/p". As 
illustrated in Fig. 2(a), P{1) is weighted towards smaller 
segments, with a maximum at P(0) = p. Thus for large 
impurity concentrations, the average over chain segments 
is dominated by contributions from the smaller clusters. 
On the other hand, in the dilute limit {p ^ 1) P{1) be- 
comes flatter, and thus larger segments contribute appre- 
ciably to the average. 

The average dynamical structure factor can conse- 
quently be calculated from 



(5) 



where the sum is taken over all contributing chain seg- 
ments with a non-vanishing Fourier component at q=7r. 
In the case of large impurity concentrations, this sum 
may be cut off beyond a certain cluster length Imax, be- 
cause the contributions of segments with I > Imax are 
exponentially suppressed by P{1). Therefore, as long 
as Imax is smaller than the maximum cluster which can 
be studied with present computer capacities (< 36 for 



Heisenberg systems), very high accuracy can be achieved 
by simply taking the sum in Eq. 5 over all available fi- 
nite systems. However, in the dilute impurity limit this 
procedure will necessarily break down, and we propose 
to use instead the fits obtained in Eqs. 1 and 2. Jl^ 

To estimate the quality of such a few-pole approxima- 
tion, we compare the spectrum S{'k, lo) obtained from 
cluster diagonalizations (Fig. 2(b)) at a rather large im- 
purity concentration (p=0.3), where Imax < 24 and this 
procedure is most reliable, with that of a three-pole ap- 
proximation using Eqs. 1 and 2 and extending the sum 
in Eq. 5 up to Z = 10000 (Fig. 2(c)). As can be seen by 
comparing the figures, the agreement between these two 
methods is very good for large p: the integrated spec- 
tral weight differs less than 0.5% at p=0.3. However, in 
the dilute impurity limit the three-pole approximation 
does account for the spectral weight due to the larger 
segments, yielding for example 14% more integrated in- 
tensity than the finite cluster average with / = 4, 24 at 
p=0.1. In the following, when discussing the evolution of 
the dynamical spin structure factor with impurity dop- 
ing, we will thus use the few-pole procedure rather than 
the average over system sizes available to exact diagonal- 
ization calculations. 
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FIG. 2. (a) Distribution function of chain segments in a 
one-dimensional system with randomly placed non-magnetic 
impurities. P{1) is shown for impurity concentrations p = 
0.1, 0.3, and 0.5. (b) Dynamical structure factor at q=7r 
for an antiferromagnetic spin-1/2 Heisenberg chain with p = 
0.3, calculated from an average over finite cluster spectra with 
1 = 4,..., 24. (c) Same as (b), but calculated using a three-pole 
approximation. 

In Fig. 3 the evolution of S{7r, lu) is shown as a function 
of p. There are three main features that occur in S^tt, uj) 
upon randomly introducing vacancies into the chain: (i) 
the integrated spectral weight decreases rapidly with in- 
creasing p, (ii) a pseudo-gap develops at small frequen- 
cies, (iii) at larger frequencies (of order J) a discrete peak 
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structure emerges, dominated by the contributions of the 
smallest cluster segments. 

In order to discuss the low-energy features (i) and (ii), 
it is useful to consider the continuum limit of Eq. 5 since 
they are dominated by the contributions of the larger 
chain segments. Wo emphasize, however, that the dis- 
creteness of the peak structure at higher energies - which 
we consider to be a major characteristic - cannot be 
accounted for in this limit. For simplicity let us con- 
sider the single-pole approximation to leading order, i.e. 
coi{l) = ai/l and Ai{l) = bi ln(l + I) (for the weight of 
the lowest pole: ai =0). Replacing the sum in Eq. 5 
with an integral it is then found that 



S{n,uj) » — o-^H^ H )exp(- 



-aip 



), (uj^J). (6) 



Hence at low energies ^(Tr,^) is exponentially sup- 
pressed, and the peak of the corresponding pseudo-gap 
feature moves towards larger frequencies with increasing 
p. The frequency sum-rule in the single-pole approxima- 
tion yields ^(Tr) = / dLoS{TT,ui) « bi exp{p)Ei{p), where 
Ei{p) is the first exponential integral. This result is in 
good agreement with the p-dcpendcncc of the frequency- 
integrated spectral weight obtained within the full three- 
pole approximation shown in Fig. 4(a). 
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FIG. 3. Evolution of S{n,iu) in the antifcrromagnctic 
spin- 1/2 Heisenberg chain as a function of impurity concentra- 
tion, calculated from a three-pole approximation. In (a)-(f) 
the respective spectra are shown for p=0.0, 0.1, 0.2, 0.3, 0.4, 
and 0.5. 

From Eq. 6 it is seen that, strictly speaking, there 
should be no spectral weight at w = 0. However, if the S- 
pcaks in Eq. 3 arc given a finite lifetime, i.e. by replacing 
them with Lorentzians of width e in order to mimic the 
experimental situation, there will be some residual spinon 
density of states 5(7r, 0). As shown in Fig. 4(b), this 
quantity is exponentially suppressed with increasing p, 
and it diverges in the dilute impurity limit. 



According to our interpretation of the pseudo-gap by 
taking the continuum limit (Eq. 6), it arises due to the 
exponentially suppressed contributions of the largest seg- 
ments (corresponding to the smallest finite-size gaps) to 
the average in Eq. 5. On the other hand, at higher ener- 
gies (to w J) the discrete nature of the smaller contribut- 
ing clusters should be visible, with the largest features 
coming from the smallest clusters. Clearly, the lowest 
pole of the smallest contributing cluster at lu = J is well 
defined in Figs. 3(b)-(f), advancing into the dominating 
feature of the spectrum for large p. At smaller impurity 
concentrations (Figs. 3 (b)-(d)) the lowest pole of the 
second smallest cluster is also well separated from the 
low-energy continuum. At larger p, however, it merges 
into the low-energy pseudo-branch. Summarizing the re- 
sults for S{tt,u!) in the presence of randomly placed static 
vacancies, an inelastic neutron scattering experiment on 
such a system should observe a marked decrease of spec- 
tral weight with increasing impurity concentration, a low- 
energy pseudo-gap feature, and - resolution permitting - 
some split-off "high-energy" peaks due to (and indicat- 
ing the size of) the smallest segments contributing to the 
average. 
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FIG. 4. (a) Dependence of the frequency-integrated dy- 
namical structure factor at wavevector tt on the concentration 
of non-magnetic impurities, (b) Residual spinon density of 
states as a function of impurity concentration. 

In real materials, impurities do not necessarily enter 
a chain in an entirely uncorrelated manner, but may 
e.g. due to a lattice potential commensurate with the 
surrounding crystal structure favor certain distances be- 
tween impurity sites. Effects of this kind can be ac- 
counted for by refining the probability distribution func- 
tion in the above analysis. E.g. in the presence of a single 
preferred segment size Iq one finds 



p(0 



p 



2 - ex.p{-plo) 



exp{-\l - lo\p), 



(7) 



which is peaked at I = Iq, and hence weighs the contri- 
bution by segments of length Iq most highly in the av- 
erage. A spectral analysis of inelastic neutron scattering 
experiments on impurity doped systems can thus indicate 
the underlying distribution function of chain segments. 
In particular, at large p the characteristic "high-energy" 
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features coming from the segments with length Iq should 
dominate the spectral weight. 

Odd-site segments do not contribute to the dynamical 
structure factor at q = tt, but they have spectral weight 
at other momentum transfers which are contained in their 
Fourier spectrum. While even-site segments have a sin- 
glet ground state associated with a gap in their excitation 
spectrum, odd-site segments have a doublet ground state, 
and the effective low-energy spin- 1/2 degree of freedom 
in these clusters leads to a single-spinon contribution to 
S{q, io) at small frequencies. This low-energy spectral 
weight from wavevectors close to tt may "spill over" , caus- 
ing a slight enhancement of the zero-frequency spinon 
density of states ^(Tr, 0). 

If there are residual inter-chain exchange couplings J± 
( I J_i| <ti \J\) in a quasi-one-dimensional magnetic mate- 
rial, there is a transition to a three-dimensional ordered 
phase at a low temperature Tn ~ |<^J-|- |l6| However, in 
the presence of an impurity-induced pseudo-gap in the 
chain spectrum, ^(Tr, 0) does not diverge anymore at low 
temperatures, and thus is decreased with increasing 
p. In the case of a vanishing residual spinon density of 
states this transition is completely suppressed, as it is 
the case for example in weakly coupled two-leg Heisen- 
berg ladders, jl^ In this sense impurities rnay have a 
stabilizing effect on one-dimensional phases. 

As a consequence of the spin-charge separation in one- 
dimensional systems away from half-filling, the spin and 
charge degrees of freedom can be treated equivalently. 
Therefore it can be expected for a segmented chain with 
(non-localized) mobile charge carriers that the results 
obtained above for S{q,u!) hold also for the dynamical 
charge response N{q,Lo). 

To conclude, we have proposed an averaging procedure 
to account for the spectral features of one-dimensional 
systems in the presence of non-magnetic impurities. At 
small doping concentrations, the contributions of the 
longer segments become increasingly important, and the 
average over their lowest poles yields a pseudo-gap fea- 
ture in the spin excitation spectrum at low frequencies. 
At larger impurity concentrations, the spectrum is dom- 
inated by discrete peaks at w ~ J contributed by the 
shortest segments. 

When static vacancies are present in a correlated one- 
dimensional system, the formation of a pseudo-gap as 
discussed above will compete or coexist with other low- 
energy features, such as gaps due to frustrating exchange 
interactions or spin-Peierls instabilities, jl^-g^ In such 
systems it then remains to decide how much of the ob- 
served discrete low-frequency peak structure is due to 
such bulk instabilities or maybe to finite segment contri- 
butions. 
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